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0\ ■ Abstract 

^! We investigate reheating of the universe by early formation of stars and quasars 

O 

\ in the hierarchical clustering scheme of cold dark matter scenario, with perturbation 

a • fluctuations normalized by the COBE data. It is found that ionizing uv flux from OB 

- , , , stars with the abundance given by the standard initial mass function is strong enough 

' to ionize the universe from z 30 to the present epoch, if 1-2% of the collapsed 

baryons go into stars. This lessens significantly the CMB anisotropics at a small 
angular scale. Reionization also increases the Jeans mass to Miuminous ~ lO^Mo for 
z < 10, which leads to a cut-off of the luminosity function of normal galaxies on a 
faint side. A strong uv flux is expected at 2; ^ 2 — 5, and the null result of the Gunn- 
Peterson test is naturally explained. Early star formation also results in production 
of heavy elements, and the observational metal abundance sets a strong constraint on 
the photon energy injection into the intergalactic space. 
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1 Introduction 



The currently most attractive scenario for cosmic structure formation assumes the 
cold dark matter dominated universe and the initial fluctuations originating from 
adiabatic Gaussian perturbations with scale-invariant spectrum (Peebles 1982, Blu- 
menthal et al.l984). The fluctuations enter non-linear regime at redshift z < 50, 
depending upon the mass of the fluctuations and the magnitude of each peak of the 
Gaussian noise, and make gravitationally bound objects. A number of N-body simu- 
lations have shown that this is indeed a successful scenario for formation of large scale 
structure (e.g., Frenk et al. 1988) and also of galaxies (White & Frenk 1992). The 
evidence in support for this scenario has been given by the COBE-DMR experiment 
for anisotropies of cosmic microwave background radiation (CMB); the data show 
that the fluctuations just after the epoch of recombination for angular scales larger 
than a few degree are close to scale-invariant spectrum with the magnitude almost 
correctly given by the CDM model of structure formation (Smoot et al. 1992). It is 
also pointed out that the CDM dominated universe predicts correctly the epoch and 
the abundance of bright quasars with the assumption that they were formed from 
high standard deviation peaks of the fluctuations (Efstathiou & Rees 1988). 

At a more quantitative level the standard CDM model does not quite flt the precise 
perturbation spectrum given by COBE-DMR and that from correlation functions of 
galaxies (Efstathiou, Bond & White 1992; Davis, Summers & Schlegel 1992), which 
motivates some authors to consider modiflcations of the model (Adams et al. 1992; 
Davis, Summers & Schlegel 1992; Klypin et al. 1992; Efstathiou, Bond & White 
1992). The original model, however, seems to be at least a good enough starting 
ground for a consideration of the history of the universe. 

In the CDM scenario structure formation is hierarchical. A small object becomes 
a gravitationally bound system flrst, and the mass of the bound systems gradually 
increases by merging. In this scenario formation of objects with a speciflc mass is 
not instantaneous, but statistical, and the abundance of bound objects with a speciflc 
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mass varies as a function of redshift. 

When bound systems are formed, it is likely that a substantial fraction of baryons 
form stars, if the cooling time is shorter than the dynamical time. The photons 
emitted from OB stars would then ionize the universe, which would affect formation 
and evolution of cosmic structure, as first pointed out by Doroshkevich, Zeldovich 
and Novikov (1967) and then discussed by Couchman and Rees (1986). In this paper 
we examine detailed thermal history of the universe in the epoch at ^ < 100 taking 
account of the effect of star formation by explicitly solving the evolution equation for 
thermal history, and discuss the consequences on formation of galaxies and background 
photons in the universe. We also take into account the uv flux from early quasars 
expected in the CDM model. 

Another interesting consequence of early star formation is that it results in pro- 
duction of heavy elements, as needed from the metal observation of globular clusters. 
Actually the observed metal abundance in population II objects sets a strong con- 
straint on the fraction of baryons which go into to stars, and hence on the energy 
injected into the intergalactic space. 

In our work we take the Press- Schcchtcr theory (Press & Schechter 1974) for 
hierarchical clustering, since it is known to describe well the clustering observed in N- 
body simulations (e.g., Efstathiou & Rees 1988). We flx the normalization of the 
Gaussian field by the data given from the COBE-DMR experiment. We confine 
ourselves to a flat universe with Q = 1 for simplicity. We assume that a certain fraction 
of collapsed baryons become stars, in so far as the cooling condition is satisfied. We 
also assume the initial mass function as observed in the Milky Way today (Shapiro & 
Teukolsky 1983), and adopt blackbody spectrum with temperature specified by the 
temperature-mass relation for population II stars. Quasar formation is normalized to 
the observed abundance at z — 2. The only free parameters of our calculation is the 
fraction of baryons which form stars apart from the Hubble constant. 

The aspects which we pay a particular attention in the present work are as follows: 
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Whether the universe is reionized and whether it recombines again at low redshift is 
of our prime interest. If the universe would be reionized early enough comparable 
to a unit Thomson optical depth, this lessens anisotropics of CMB at an angular 
scale smaller than a degree, which would make them consistent with the limit from 
the South Pole II experiment (Gaier et al. 1992). Ionization of the universe would 
increase the Jeans mass, which otherwise is very small compared with the galaxy mass 
today, to the order of irregular galaxies. The exact Jeans mass limit depends on how 
efficient is ionization. The evolution of residual photons is also of significant interest. 
For z < b the Gunn-Peterson test (Gunn & Peterson 1965; Steidel & Sargent 1986; 
Jenkins & Ostriker 1991; Webb et al. 1992) requires that the universe be ionized to 
a high degree. This means that a reasonable amount of the uv flux must exist all the 
time for the relevant redshift range to keep the universe ionized. The necessity of uv 
flux at around 2; ~ 4 — 2 is also motivated by the presence of Lyman a clouds; the uv 
flux heats up the clouds so that they do not collapse (Sargent et al. 1980; Ostriker & 
Ikeuchi 1983; Rees 1986). 

The effect of reheating has been studied by a number of authors (Couchman 1985, 
Stcbbins & Silk 1986; Couchman & Rees 1986, Bartlett & Stebbins 1991; Barcons et 
al. 1991; Fukugita & Kawasaki 1990; 1993; Gnedin & Ostriker 1992; see Cen et al 
1990; Cen & Ostriker 1992 for a different class of the reheating model) in different 
contexts and models. In particular Couchman & Rees (1986) studied quahtatively the 
possibility of reheating in the hierarchical clustering scenario. We shall study more 
quantitatively the thermal history of the universe, with the strength of the fluctuations 
determined from the CMB anisotropy observations. Reheating in the hierarchical 
clustering was also considered by Sasaki et al. (1993) However, these authors take 
account only of the energy liberated from the gravitationally binding objects. Such 
scenario basically does not ionize the CDM universe. The energy injection from stars 
is crucial for reionization. 

In sect. 2 the basic model of the present calculation is described. The method 



4 



of the calculation for the thermal history is given in sect. 3; the thermal processes 
that we considered are described in detail. In sect. 4 we present the results of our 
calculation and discuss the meaning and implications of the results. Sect. 5 is given 
to a conclusion. 

2 The model 

2.1 Press-Schechter theory 

In the Press-Schechter theory the comoving number density of non-linear objects 
within the mass range between M and M + dM is given by 

where po is the mean comoving mass density, 5c is the overdensity threshold for the 
collapse (= 1.68 corresponding to the prediction of spherical collapse model), Di{Z) 
is the Peebles function for the growth of the perturbations {Di[z) = (1 -|- z)~^ for 
0=1), and a'^{M) is the Gaussian variance given by 

a\M) = —j^ P{k)eM-rlk')edk, (2) 

with rM — (M/po)^^^/'\/27r the width of the Gaussian filter and P{k) the power 
spectrum of density fluctuations. We use the fitting formula for the CDM power- 
spectrum given by Bardeen, Bond, Kaiser & Szalay (1986): 

P(k) = A\T(k)\% (3) 

_ ln(l + 2.34g) 
^^"^ - 2:34^ 

X [1 + 3.89g + (16.1g)2 + (5.46?)^ + (6.71g)^]-^/^ (4) 
k 

where h is the Hubble constant Hq in units of lOOkmsec"^ Mpc"^). This gives a 
spectrum quite similar to that given in Efstathiou, Bond & White (1992), but with a 



slightly more power for small scales. The normalization A of the spectrum is fixed by 
the COBE-DMR data using the relation (Efstathiou, Bond & White 1992): 



Qrms / 5 \^^^ [ Hq 
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(6) 



where Tq is the present temperature of CBR and Qrms is the fluctuation amplitude at 
~ 10° (Smoot et al. 1992): 

grms = 16.7 ± 4.6 //K. (7) 

This normalization fixes N{M, z) now uniquely. 

The fraction OcoU of the total mass of the universe which collapses into the bound 
object with mass greater than M is shown in Fig. 1(a). The effect of reionization is 
not taken into account here. The Jeans mass before reionization is given by (Bond & 
Szalay 1983) 

Mb,j = 1.4 X lQ^MQVtBVL-^''^h-^{l + 100/^)-^/^ (8) 

which stays at a constant value 1.4 x 10^Q~^/^Mq for z > 100 due to residual ion- 
ization left over the recombination epoch. The total mass corresponding to this Jeans 
mass is 3 X IO^Q'^^'^Mq (or less for a smaller z). 

Another necessary condition for continuing collapse is that the cooling time is 
shorter than the dynamical time (Blumenthal et al. 1984), 

/ V ^ -0.625 s -1.04 /I I -2.75 

Ms,>3.7xio^nhr'-'''^] m r-^] (9) 



.lo-v \n J V 10 

where Ye ~ IQ-^n^^n^/^h-^. At z = 50 the hmiting mass is M ~ lO^M© and 
decreases as a function of redshift. Therefore we conclude that the first objects which 
collapse around 2; 50 should have a mass of lO^M©. 

In the following calculation we take account of the effect of reionization on the 
Jeans mass, which is written 

Mb,j = 1.4 X 10^MoQsQ-=^/2/i-i i-^j , (10) 
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where is the electron temperature, is the temperature of the cosmic background 
radiation and /x is the mean molecular weight. When the universe is reheated, 
becomes much larger than and increases. Therefore, Jeans mass Mb,j increases 
whereas Mb,c decreases. We take the object with m.&yi[{Q,/Q,B)MB,c, {^/^b)Mb,j\ < 
M < IO^^Mq as being collapsed. For M > lO^^M© the cooling time becomes longer 
than the Hubble time again, and galaxies are not formed (Blumenthal et al. 1984). 

2.2 Early formation of stars 

When bound objects are formed, we assume that a constant fraction / of the baryonic 
gas goes into stars in so far as mass of the objects is greater than the Jeans mass and 
smaller than the maximal mass of galaxies (= IO^^Mq), and the cooling time is shorter 
than the dynamical time. We set / = otherwise. We take the baryon density to be 
^B — 0.05/i~^ and Qb/^ to be universal. We also assume that the mass spectrum of 
stars is proportional to the present-day initial mass function (IMF), ^imf, for which we 
adopt the formula given by Shapiro & Teukolsky (1983). We take the temperature- 
mass relation for population II stars (e.g.. Bond, Carr & Hogan 1986): 



Stars with mass smaller than IOMq produce very few ionizing photons. Therefore it is 
sufficient if we only consider the uv flux from stars whose mass is greater than IOMq. 
A star with mass Mg produces the radiation energy SgMg in the main-sequence time 
^MS) where eg and ^ms are given by 



where X = 0.76 is the hydrogen mass fraction. Since the main-sequence time of 
massive stars with Mg > IOMq is shorter than the cosmic time for z < 50, we may 
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(11) 




(12) 



and 




(13) 



assume that the uv photons are produced instantaneously after the formation of the 
stars. The production rate of uv photons is then 

C^L ^ (^) ■ (") 

where B{e^,Ts) is the blackbody spectrum normahzed to / de^B{e^) = 1. 

Stars with mass greater than AMq eject heavy elements. Following Carr, Bond & 
Arnett (1984) we assume that the fraction of mass ejected as heavy elements is 

Zei = 0.5 - for 15 < M/Me < 100 (15) 

Zej = 0.1 for 8 < M/Mq < 15 (16) 

Zej = 0.2 for 4 < M/Mq < 8. (17) 

Then evolution of the metallicity Z is given by 




In our work we take / = 0.02 which, as we find, is about the upper limit from the 
condition that metaUicity of population II stars does not exceed 10~^. 

We show in Fig.l (b) the mass fraction of stars both with and without taking the 
effect of reionization into account. The figure indicates that a half of the stars are 
formed before z — 10. This is because we assumed that stars are made always at a 
constant fraction of collapsed baryons with the cut-off at M < IO^^Mq. In reality 
we should have more stars formed stationally from the collapsed gas, and expect that 
our approximation would not be good at low rcdshift, after galaxies arc formed to 
their present form. Fig. 1(c) shows the evolution of the emitted radiation energy from 
stars. In principle, our scheme has a double counting problem in that the baryons 
which have gone into stars when a small object collapses might be used again to make 
another stars when a larger object is formed. Since the fraction which goes into stars 
are very small (~ 10~^), however, the double counting effect is actually negligible. 
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2.3 Formation of quasars 

We take account of the effect of quasars following Efstathiou & Rees (1988). Using 
the Press-Schechter theory, we estimate the comoving number density of luminous 
quasars as 

A^Cuasar(^) = f [ ( ^^^^''A dMdt (19) 

Jmax[t-tQ,0] JMmin \ Cft J 

where tq is the lifetime of quasars which is assumed to be much shorter than cosmic 
time t, and M^am is the minimum mass of quasars. Wc take Lq^min = 10'^'' erg scc~^ as 
a cut-off luminosity of quasars and estimate Mmin from the mass-luminosity relation: 

MiL,> ^ 2 X lO^'M, (^) (^) , (21) 

where Fq is the fraction of the collapsed matter which becomes quasars, and eg is 
the efficiency defined as the fraction of the rest mass energy that is converted into 
radiation. The quasar number density depends on and cqFq. Wc take tq = 10^ 
year and fix eqFq to be 10~^ to give the observed comoving number density of quasars 
brighter than Lq = 2.5 x lO^^/i^^ergsec""^ at z = 2, 

A^qua^ar ^ 1-5 X IQ-^/i-^MpC-^. (22) 

The uv photon production rate is given by 

\ / quasar ' \ / 

where -S'Q(e-y) is the spectrum of the emitted photons which we take to obey a power- 
law of the form (Miralda-Escude & Ostriker 1990; model QS2) 

(e^ < 10.2eV) 

oc <; , (24) 



er^-^ (e^>10.2eV) 



We normalize as / de^Sq = 1. 



The abundance of quasars are given in Fig.l (d). The first quasars are formed 
around z ~ 7 — 8. The comoving number density reaches maximum at z ~ 2.5 and 
then it dechnes. The radiation energy from quasars is shown in Fig. 1(c). 

2.4 Heating by collapsing bound objects 

The possibihty was discussed that hberation of kinetic energy by radiation coohng 
during the collapse of galaxies might ionize the universe (Hogan 1980). If free- free 
cooling or recombination cooling is the dominant cooling process, the produced pho- 
tons might have energy high enough to ionize the intergalactic medium. This possi- 
bility was reconsidered within the context of hierarchical clustering theory recently 
by Sasaki et al. (1993). We can show, however, that this scenario does not work for 
the following reasons: for this scenario to work, free-free or recombination cooling (we 
denote the cooling times as Tff and Tree) must be faster than Compton cooling (tq), 
expansion cooling {ju — 1/H) and line cooling (rune) that follows coUisional excita- 
tion. The conditions {ts or Tree) < ("^c and th) are satisfied only for low temperature, 
T < 2.5 X lO^K, as seen in Fig. 2 where the critical temperature for tc = etc. are 
plotted as a function of redshift. On the other hand, we can show that, at such a low 
temperature, line cooling dominates over recombination and free- free cooling (Fig. 2). 
This means that liberation of kinetic energy by radiation cooling in binding objects 
is not efficient enough and does not ionize the universe. We give this argument more 
quantitatively in Appendix A. 

The situation gets worse if we include line cooling due to helium. On the other 
hand, the constraint might be weakened, if the inhomogeneity effect of baryon den- 
sity in collapsed objects that enhance recombination (free-free) cooling is taken into 
account, as discussed by Sasaki et al. (1993). In any case, however, only a very 
narrow range of temperature may be allowed to produce ionizing photons, even if 
this enhancement is included. The contributions from massive stars and quasars are 
far more important in reionizing the universe. We conclude that we can neglect the 
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production of ionizing uv photons from collapsing objects. 



3 Thermal history 
3.1 Method of calculation 

The evolution equations are solved numerically for the ionized fraction of hydrogen 
and helium atoms n(HII)/nH, n(HeII)/nHe and n(HeIIl)/nHc, the photon spectrum 
n^(e^,t), and the electron temperature Tg with all relevant thermal processes taken 
into account. The integration is made with the implicit Euler method. The detailed 
expressions of thermal processes used in our work are given in order in the following 
subsections. 

We determine the spectral distortion of CMB by solving directly the Kompaneets 
equation (e.g., Fukugita & Kawasaki 1990). Since reheating occurs at low z and 
the temperature is not high in the present model, the distortion is described by the 
Zeldovich-Sunyaev spectrum and its magnitude is characterized well (within the ac- 
curacy in solving numerically the Kompaneets equation) in terms of the Compton 
t/c-parameter: 



3.2 Ionization of H and He 

The abundances of HI, HIT, Hel, Hell, and Helll are determined by the balance 
between the ionizing processes and the recombination processes. We treat a hydrogen 
atom as a two level system {1S,2S + 2P) plus continuum following Peebles(1968), 
Matsuda et al. (1971) and Jones & Wyse (1985); evolution of the HII fraction is 
determined by 




(25) 




+ 



i?2cn(Hl) K{a2,minl + y\ri(Hl)e-'^'^«/'=^^) + g-^'^-/'^^^ 
^ 1 + K{a2,minl + i?2cn(HI) + An(Hl)) 



(26) 
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where nn = n(HI) +n(HII), A(= 8.227 sec~^) is the two-photon decay rate from 25, z/q, 
is the Lyman-o; frequency, K — c^/(87ri/^)(a/a), q;2,hii is the recombination coefficient 
to 25" + 2P level and i?2c (Ric) is the ionization coefficient from 25" + 2P (IS). i?2c 
and it!ic are given by 

R2c = l2c + P2,mne+ / ca2{,mn^{ej)de^, (27) 

-Ric = /3i,Hi?T'e+ / cai{^mn^ie^)de^, (28) 

where /3i,hi {P2,m) is the coefficient for collisional ionization from IS* (25* + 2P) level, 
72f = ci;2,Hii(27)(2meA;r^)^/^e~^'^®^/^i' is the photo- ionization coefficient due to back- 
ground photons, (Tif^Hi (c2f,m) is the photo-ionization cross section from 15" (2S'-|-2P) 
and n-y(e-y) is the spectrum of uv photons; we take ei^ni = 13.6eV, e2,Hi = 10.2eV. 
Explicit expressions for the coefficients used in this work is summarized in Appendix 
B. 

Since helium is a minor component, we may treat it as a one- level system. The 
time evolution of HI, HII and HIII is then given by 

d /n(Hell)\ n(Hel) 



dt \ nue 



de^aYyi^Ueicn-fi^-/) 



n(Hell) n(Helll) 

- aueline \- OiHelune 

^He ^He 

^(Hell) 

- Cneiirie- (29) 

^He 



d /n(HeIII \ n(HeII /• , 

— = / de^crbf,Heiicn^(e7) 

dt \ riHe J riHe JeHeii 

n(Hell) n(Helll) 

+ PHeime OiHeUine , (30) 

riKe "-He 

where nuc = ^^(Hel) + n(Hell) + n(Helll), P, a and ^ are the coefficients for the 
collisional ionization, recombination, and dielectronic recombination, respectively; cxbf 
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is the photoelectric ionization cross section, enei (cHeii) is the ionization energy of Hel 
(Hell) (see Appendix B). 



3.3 uv spectrum 

The time evolution of the photon spectrum is determined by solving 



dt ^ \ ^^^7 



1 . . . , 

+ — Uff + Jfb,HII + Jfb,HeII + Jfb,HeIIlj 



— (cif,HI + C2f,HI + Cbf,HeI + Cbf,HeIl)c% 

(^) ' 

\ / star \ / quasar 

where js is the emissivity of the free-free process, jVb.Hii, Jfb.Heii and jfb,Heiii are that 
of the free-bound process due to HII, Hell and HeHI, respectively. In the present 
consideration the electron temperature is relatively low ( < lO^K), and hence the uv 
photon production due to the free-free process is not important. 

The photons produced by recombination process have a complicated spectrum. 
We treat HI and Hell as two level systems {IS, 2S + 2P) and Hel as a three level 
system (1^5*, 2^5*, 2^P). When a free electron is captured by HII, Hell or Helll, the 
emmisivity is given by 

Jfb,ij = -^aijn{i)ne exp(-(e^ - eij)/kTe), (32) 

where i =HII, Hell and Helll, and j represents the energy level; aij is the recom- 
bination coefficient and e^j is the binding energy of the j-level of an i-atom. Unless 
an electron recombines directly into the IS level, a monochromatic (line) photon or 
two photons are subsequently emitted from an excited atom at every recombination 
process. The photon emissivity is 

ifb,i,i,ime = e^aijn{i)neS[e^ - (e^ - e^j)], (33) 
13 



for a monochromatic photon emission, and 

jfb,i,i,two = 2e^aijn{i)nJ[e^ - (e^ - eij)/2], (34) 

for the two-photon process.^] 

In case of dielectric recombination of Hell, two electrons are in excited states after 
recombination, and they go down to the ground state by emitting photons. We treat 
this process by assuming that the first electron goes to the ground state emitting a 
photon with energy corresponding to (25" + 2P) level of Helll and the second electron 
emits a photon with energy corresponding to 2^P level of Hell. Therefore, 

Jfb,Heii,d = e^eHeiin(Hen)ne[(5(e^ - 40.81eV) + (5(e^ - 21.2eV)]. (35) 

3.4 Electron temperature 

Evolution of the electron temperature is given by 
3 d f kTenB\ 



2dt\ /i 



n{i)c i (e^ - ti)n^(Tu,ide^ 



i=HI,HeI,HeII 
i=HI,HeI,HeII 
j=HII,HeII,HeIII 

- u;Heii^e^(HeIII) 

- Y i'inenii) 

i=HI,HeI,HeII 

- Ac 

- Os [n(HII) + n(HeH) + 4n(Heni)]ne 
15 d / kT^nB\ 



(36) 



2 a \ fi J 

where Q is the collisional-ionization cooling coefficient due to atomic state i, rji is 

the recombination cooling coefficient, c^Heii is the dielectronic recombination cooling 
^Here we neglect the effect of the energy distribution for the two-photon process. 
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coefficient due to Hell, ipi is the collisional excitation cooling coefficient due to i, Ac 
is the Compton coohng rate and 6s is the free-free coohng coefficient (see Appendix). 
The last term of (^) represents expansion cooling. In the actual calculation we also 
include Compton heating by scattering of high energy photons off electrons, although 
it is not important since the reheating occurs at low redshifts {z < 100) and the 
energy of photons is small ( < keV) in our model. Compton cooling, among other 
cooling processes, is most important at z > 5; the electron temperature is basically 
determined by the balance between Compton cooling and photoelectric ionization 
heating. At a lower redshift expansion cooling dominates. 

4 Results and implications 

The results of our calculation are presented in a few panels of Fig. 3. Fig. 3(a) shows 
a fraction of HI and HII; the universe is reionized at 2 ~ 30, shortly after the ffist 
objects with mass of ~ IO^Mq collapse. By this epoch the relative abundance of 
the baryons collapsed into stars is Qstar/^B — 10~^ and the photon energy density is 

^ 10^^'^. Since then ionization is kept at a high degree; the HI fraction decreases 
below 10~^ by 2; ~ 20 and stays at a value < 10^'' to the present epoch, which explains 
the very strong limit on neutral hydrogen from the Gunn-Peterson test (Gunn & 
Peterson 1965; Steidel and Sargent 1989; Jenkins & Ostriker 1991; Webb et al 1992). 
If we would switch off quasars, the hydrogen tends to recombine towards 2; 1 — 2 
[n(HI)/nH ~ 10~^ at z = 0]. However, the uv flux from quasars, which appears as a 
kink in the curve of HII around z ^ 6 when they turn on, ensures continuous high 
ionization to z = 0. We remark that the absorption due to possible Lyman a clouds at 
high z might reduce the ionizing flux and lead to more HI (Miralda-Escude & Ostriker 
1992; Madau 1993), which we ignore in the present work however. 

Whether anisotropies with angular scale smaller than a few degrees may be less- 
ened in our reheating scenario is an interesting issue of the model. The Thomson 
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optical depth calculated as 

poo 

tt — dtca^rif. (37) 

^0 

is Tt — 0.22 for the ionization history shown in Fig. 3(a) (the effect of helium is also 
taken into account). While this optical depth is not so large, a calculation made 
by Sugiyama, Vittorio & Silk (1993) indicates that the Doppler peak of the Fourier 
component of CMB correlation function q at £ ~ 200 decreases by a factor 2 and also 
the power for ~ 10' by a factor of 3 for reionization z ^ 30. Their calculations 
also show that anisotropics at 1 degree scale are appreciably smaller than 1 x 10~^, 
and are perfectly allowed by the limit from the South Pole 11 experiment (Gaier et al. 
1992). 

The Compton ijc parameter is also plotted in the same figure (Fig.3(a)). The 
final value is 2 x 10~^. This small yc is understood by the fact that ionization takes 
place at low redshift with low energy photons, which keep electrons to stay at a 
low temperature; hence only a small fraction of energy goes into CMB by Compton 
cooling. This contrasts the model of Cnedin & Ostriker (1991), where high energy 
photons are injected from accretion activity of massive black holes which are formed 
at a very high redshift; in their model yc is as high as 10~^, though the predicted value 
is perhaps flexible. We should also add a remark on the difference of our model from 
the model with heating by shock waves; the latter heats up electrons to 1-lOkeV and 
leads to i/c-parameter as large as > 10^^ (Yoshioka & Ikeuchi 1987); shock heating is 
allowed only at small z ( for this class of model, see Cen & Ostriker 1992). 

The metal abundance Z (Fig.3(b)) closely follows the star abundance. For z < 10, 
Z stays at (0.4 — 1) x 10~^, which is a typical metal abundance of population II stars. 
This means that our fractional star-formation parameter / ~ 2 x 10~^ is close to the 
limit; namely, the observed metal abundance limits the amount of energy injection 
from stars. 

Ionization history of helium atoms is exhibited in Fig. 3(c). He I is ionized to He II 
{AE — 24.6eV) at the epoch of hydrogen reionization, and n(HeII)/nHe — 10~^ after 
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this redshift. He II is then fully ionized to He HI [AE = 54.4eV) when quasars turn 
on. Since the starlight can not ionize He II, the abundance of He HI is taken as an 
indicator for the importance of the quasar hght (Mo, Miralda-Escude & Rees 1993). 

In Fig. 3(d) we plot the electron temperature, the baryon Jeans mass (10), and the 
critical mass for the cooling condition (9). We note that the electron temperature stays 
between 3 x 10^ — 2 x lO^K for all period of z < 30. The fact that the universe is fully 
ionized despite this low temperature can be understood by the slow recombination 
rate. 

One of the most conspicuous effect in the reionized universe is a large Jeans mass at 
low redshift (Couchman & Rees 1986). It gradually increases from 10^ Mq for z < 100 
to 1 X IO^Mq for z < 6. This means that normal galaxies formed at 2; < 10 should have 
a mass function with a peak around Miuminous ~ 10^ — lO^'', or a luminosity function 
with a peak at a S magnitude Mb ~ — 19mag rather than obey a Schechter function. 
Observationally, little is known about the faint end of the luminosity function of field 
galaxies. The only information available on the faint end is from galaxies in the Virgo 
cluster; Binggeli, Sandage & Tammann (1985) obtained a complete galaxy sample 
down to B=18mag, and they have shown that the luminosity function of spiral and 
elliptical galaxies are of the Gaussian type with peaks at Mb — — IS.lmag (dispersion 
a — l.Smag) and Mb — — 18.2mag {a — l.Tmag), respectively (Sandage, Binggeli & 
Tammann 1985). [Here we used ((m — M)o) ~ 31.4 for spiral galaxies and 31.0 for 
elliptical galaxies; see the discussion in Fukugita, Okamura & Yasuda (1993).] The 
Jeans mass curve shown in Fig.3(d) suggests that very small galaxies, such as dwarf 
spheroidals, and globular clusters are very early objects which collapsed at 2; > 10 
and have survived merging to larger systems. 

The critical cooling mass is smaller than the Jeans mass in the reheated universe, 
in contrast to the case without reionization, for which generally Mb,c > Mb,j- 

The resulting photon spectrum is given in Fig.4 for z = 5,3 and 0. The relative 
contributions from stars and quasars are obvious in this figure. We show in Fig. 5 
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the strength of ionizing flux deflned by 



•'El 



em 47rC 



~ de. 




Clf,Hl(eHl) 



(38) 



as a function of redshift. The ionizing flux rises very sharply at z ~ 30, the epoch 
of reionization. As we noted earher the flux at low redshift may have substantial 
uncertainties, arising from our underestimate of stationally star formation activity in 
normal galaxies and from the neglect of the absorption by Lyman a clouds. Never- 
theless, the value of J around z — 5 — 3 {J — {3 — 5) x 10~^^erg cm~^ s~^ sr~^ Hz~^) 
is interesting, since it is about the value required to keep Lyman a clouds from their 
collapse (Sargent et al 1980; Ostriker & Ikeuchi 1983; Rees 1986; Bajtlik, Duncan & 
Ostriker 1988). 

5 Conclusion and discussion 

We have shown that OB stars which would have been formed shortly after the collapse 
of very early objects, as expected in the standard CDM model, ionize the universe 
completely at the epoch as early as z 30. This reionization signiflcantly lessens 
the anisotropies of the CMB at an angular scale smaller than 1° and removes the 
marginal conflict between the COBE-DMR data and the limit from the South Pole II 
experiment in the CDM model (e.g., Gorski, Stompor & Juszkiewicz 1992), without 
invoking the aid of the tensor perturbations (Crittenden, R. et al. 1993; Davis, R.L. 
et al. 1992; Dolgov & Silk 1992). Reionization has also an interesting imphcation, 
that luminosity function of normal galaxies formed at z < 10 should have a cutoff on 
the faint side, which roughly corresponds to the luminosity of dwarf galaxies. Smaller 
bound systems such as globular clusters and dwarf spheroidals must be the objects 
formed very early in the universe. 

Another advantage of the present scenario is that the metal abundance of globular 
clusters can be explained naturally. The observed metal abundance sets an upper 
limit on the photon energy that is emitted from stars and injected into the universe. 



18 



Wc have shown that the null result of the Gunn-Peterson(GP) test is satisfied by 
reheating in contrast to the case of shock wave reheating as advocated by Cen and 
Ostriker 1992; with shock heating alone ionization is not so efficient and the GP test 
is not satisfied by many orders of magnitude. The GP test is not satisfied in a heating 
model with very early formation of massive black holes (Gnedin & Ostriker 1992); 
the light emitted from an early epoch is redshifted away and has no ionizing power at 
low redshift. In order to keep hydrogen highly ionized, the ionizing uv fiux must exist 
during the relevant epochs. This is satisfied in our model by successive formation of 
bound objects. We also expect the ionizing fiux as strong as is required to confine 
Lyman a clouds by continuous ionization. 

In our reheating scenario the Compton ijc parameter takes a very small value, 
much smaller than could be detected by the COBE-FIRAS experiment (Mather et al. 
1993). 

The most interesting test for the present case is the CMB anisotropics at a small 
angular scale. If anisotropics would be observed, for instance, at a level of AT/T 
1 X 10~^ reheating must not have happened before z fa 10 — 15, which means that 
the baryon fraction that goes into stars should be smaller than / ~ 10~^ in our 
terminology. 

All of our calculation, except for the normalization of the fiuctuations, scales as 
Vlh? (and Qs/i^), and the normalization of the fiuctuations depends very weakly on 
Q (f2~°-^). Therefore, the reheating scenario changes little if high value of the Hubble 
constant is taken in a low density universe. If Vlh^ <^ 1, however, the epoch of 
reheating is considerably delayed. 

Our calculation presented here is certainly far from complete for thermal history 
and evolution of galaxies after they are formed. We have to carry out a calculation 
with a more astronomical treatment to take proper account of the effects. 
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A Quantitative argument for heating by collapsing 
bound objects 



When the matter collapses, baryon gas and dark matter particles are virialized and 
distributed isothermally with radius Tyir and temperature Tyir, which are given by 



M 



nl/3 



= 1.7 X lO-^Mpc ( (l + z)-^(Q/i2)-i/3^ 



(39) 



and 



_ GrUpMii 
'''' ~ 3r • 

Ol vir 

= 1.04 X 10-^K /X j (l + z)(Q/i^)V3, (40) 

where Pvir is the mean mass density of the virialized object and we take Pvir = 180p 
from the spherical collapse model. The temperature of the baryon gas decreases by 
radiative cooling processes, which convert the kinetic energy of baryon gas into ra- 
diation. If free-free cooling or recombination cooling dominates over other cooling 
processes, the produced photons could have energy high enough to ionize the inter- 
galactic medium. 

This scenario requires that free- free or recombination cooling must be more efficient 
than Compton cooling, expansion cooling and line cooling. For fully ionized plasma, 
the cooling times for the free-free and recombination processes are 

1.66 X 10^^ sec {T/KY'^{1 + z)-\nBh'')-\ (41) 



Tflf = 



and 
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T 

~ 1.19 X 10" sec i- 

K 



0.7 



T 



106K 



0.7- 



2\-l 



(42) 



where the coohng rates and Arec are given by 



As = 1.42 X 10'^V(^e/K)^/^(n(HII) + 4n(HeIII))ne erg cm"^ sec"\ (43) 



. 6.5 X 10- (I)"' (^) 

xn(HII)ne erg cm~^ sec~^, 



-0.2 



1 + 



0.71 -1 



(44) 



with He, n(HII) and n(Helll) the number densities of electron, HII and Helll, and 
Qs the gaunt factor of ~ 0(1). That the coohng time or Tree be shorter than the 
Compton coohng time tq = 3.9 x 10^^(1 + z)~^ sec leads to the conditions: 



Tc 

Tree 
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65 



/ 7^ \-i/2 



lO^K 

T 
lO^K 



0.0125 



-0.7 



> 1, 



(45) 



0.0125, 



> 1. 



(46) 



The cooling times or Tree should also be shorter than the dynamical time th 
1/ H = {a/a)~^ to make the cooling processes effective. 



H 



-1 



vio^k; 



(1 + 



2\ -1/2 



0.0125/ V0.25 



> 1, 



(47) 



Tr, 



O.T 



(1 + ^) 



3/2 



-1/2 



> 1. 



(48) 



V0.0125 ; V0.25 / 

These condition (|45|)-(|48|) are satisfied simultaneously for temperature T < 3 x lO^K, 
as seen in Fig. 2, where the critical temperature for tq = Tg etc. are plotted as a 
function of redshift. 

Another condition to produce efficiently ionizing uv photons is that the recom- 
bination cooling rate must dominate over the rate of line cooling of hydrogen atoms 
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that are excited coUisionally. The latter rate is given by 

^/^l"^ / 1.18x10^^ 



Aline = 7.5 X 10 
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exp 



nen(HI) erg cm ^sec ^, 

(49) 



so that the ratio Aune/Arec is 

Aline 



A. 



2.9 X 10^ 



-0.3 



1 + 



1/2N 



X 



V106K/ 



o.7^ 



n(HI) 



exp 



-1.18 X lO^K^ 



(50) 



The condition Aune/Arec < 1 means that the fraction of HI be small enough. Fig.6 
shows the contour of Aune/Arec = 1 in the plane of T versus [n(HI)/n(HII)]; the 
production of photons requires that the parameters should lie below this contour. On 
the other hand, n(HI)/n(HII) is determined by ionization equilibrium, as 



n(HI) 
n(HII) 



X 



4.95 



1 + 



-1.2 



1 + 



lO^K 



0.5 



0.7- 



exp 



1.58 X lO^K' 



(51) 



Vio^k; 

The curve showing the equilibrium condition is also plotted in Fig.6. One can see that 
the temperature should be greater than 3x lO^K for recombination cooling to dominate 
over line cooling. However, this is just opposite to the condition for recombination 
cooling being faster than Compton cooling. This proves that liberation of energy in 
bound objects does not efficiently produce ionizing photons. 

B Ionization and recombination coefficients 

The ionization and recombination coefficients used in our calculation is summarized 
in this Appendix. The coefficients and cross sections used in our paper are taken from 
Black 1981; Matsuda et al. 1971; Menzel & Pekeris 1935; Spitzer 1978. We adopt 
the high temperature correction given by Cen 1992 for the coefficients of processes 
involving a free electron and an orbital electron. 
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B.l Collisional ionization 



(a) HI (n = 1 ^ free): 



Pi,m = 5.85 X 10-^^r^/^(l + (T/10^)°-^)-^ exp(-1.578 x W^T) cm=^sec-\ (52) 

(b) HI (n = 2 ^ free): 

P2,m = 5.46 X IQ-^^T^/^ exp{-3.882 x 10^ /T) 

X [19.98 - 5.89 X lO^^T - 2.81 x lO^T'^ + BAA x lO'^T''^] cm^ sec-\(53) 

(c) Hel: 

/^Hei = 2.38 X 10-^^r^/2(l + (r/10^)°-^)-^exp(-2.853 x lOV^) cm3sec-\ (54) 

(d) Hell: 

Pneii = 5.68 x 10-^^X^/2(1 + (T/10^)°-^)-^ exp(-6.315 x lOV^) cm^ sec"^ . (55) 

B.2 Recombination 

(a) HH (free > 1): 



ami = 6.28 x IQ-'^^T-'^/^ 



T 
103 



-0.2 



1 + 



T 
105 



0.7- 



cm^sec ^, 



(56) 



(b) HII (free ^ n > 2) [Peebles 1993]: 



s -0.8 

a2,Hii - 2.6 X 10-^3 ( j^j cm3 sec'S 



(57) 



(c) HeU (free > 1): 



ttHeii = 1-50 X io-iot-°-^3^3 cm3sec-\ 



(58) 



(d) HcII (free ^n^[n>2]) 



Q;Heii23 = 9.94 X io-iir-°-^^^^ cm^sec-\ 



(59) 
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(e) Hell (free ^ l^S) 



aneiiiis = 1-32 x iq-^^T'^-^^^ cm3sec-\ 



(f) Hell (dielectronic recombination): 



^Heii = 1-9 X 10-=^T-i-5exp(-4.7 x IOVT) 

x[l + 0.3exp(-9.4 x lOV^)] cm^sec"\ 



(g) HeHI (free > 1): 

aneiii = 3.36 x IQ-'^^T-'/^ (^^^ 



T 



4 x 106 



0.71 -1 



cm 



(h) HeHI (free > 2): 



ttHeIII,2 



Q;Hem(l.ll - 0.044 In T) T < 2.2 x 10^ 
Q;Hem(1.43 - 0.076 In T) T > 2.2 x 10^ 



B.3 Photoionization cross sections 



(a) HI (n = 1 ^ free): 



(7bf,Hi = 1.18 X 10-^^e-^- 



-4(arctan 2i)/zi 



1 — g-27r/zi 



cm 



where Zi — [e^/eni — 1]^^'^. 
(b) HI (n = 2 ^ free): 



Cbf,HI,2 



1o , 3 + 4^2 5 + 4^2 e-4(arctan(22))/2 

= 1.08 X 10-i^er'^^— ^V?V^7T^ ^T^ 



where z — [e^/eni — 1/4] 2. 
(c) Hel: 



(7bf,Hei = 1-13 X 10 



1 9.775 



(d) Hell: 



10-4' 



(J6/,Heii = 7.55 X 10-^^e^ 



where Z2 = [e^/eneii - 1]^/^- 



— 4(arctan Z2)/z2 

\ — g-27r/22 



cm 



25 



B.4 Cooling Rates 

B.4.1 Collisional ionization cooling 



(a) HI: 

Chi = 1.27 x IQ-^^T^/^ 

(b) Hel: 
Cnei = 9-38 X 10-22r^/2 

(c) Hel(235): 

CHei(235)n(HeI) = 5.01 X 10-^^r-°-^6«^ 



T \ 1/2- 



T \ 1/2- 



exp(-1.58 X 10^/T) erg cm^sec" 



exp(-2.85 X 10^/T) erg cm^sec 



1 + 



X \l/2- 



105 



X exp(-5.53 X 10 /T)nen(Hell) erg sec . 



(d) Hell: 
Cneii = 4.95 X 10-22^^/2 



T \ 1/2 



exp(-6.31 X lOV?") erg cm^sec 



B.4. 2 Recombination cooling 



fa) HH: 



?7hii = 6.50 X 10' 



-27J.1/2 f 



T 
Vl03 



-0.2 



0.7- 



Vloe; 



erg cm^sec ^, 



(b) Hell: 



r^Heii = 1-55 X 10-26^0.3647 ^^.g cm-=^ sec^S 

(c) HeH (dielectronic recombination) 

u;Heii = 1-24 X 10-^^r-^-^exp(-4.7 x IOVT) 

x[l + 0.3exp(-9.4 x lOV^)] erg cm2sec-\ 

(d) Helll: 



77Heiii = 3.48 X 10-26rV2 (^)'° ' 



1 + 



4 X 106 



0.7- 



erg cm^ sec ^ 
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B.4.3 Collisional excitation cooling 



(a) HI: 

^Hi = 7.5 X 10-^^ 

(b) Hel 



1 + 



rjn ^ 1/2- 

105 



exp(-1.18 X lOVT) erg cm3sec-\ (76) 



V'Hein(Hel) = 9.10 X 10-'^T-°-^^«^ 

X exp(-1.31 X lOVT)nen(Hell) erg sec , 



^+'10^ 



(77) 



(c) Hell: 

^Heii = 5.54xlO-^^T-°-='9^ 
B.4.4 Free- free cooling 

B.4.5 Compton cooling 



rp . 1/2- 



exp(-4.73 X IO7T) erg cm'' sec"' . (78) 



= 1.42 x-27^ffTV2. 



TT^ (kT\ ( kT \ 



(79) 



(80) 
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Figure Captions 

Fig.l (a) Evolution of the density of collapsed objects with mass M greater than 
prescribed values (shown in units of Mq) for Vt = 1, h = 0.5. 

(b) Evolution of the mass fraction of stars for Ob = 0.05 and / = 0.02 {VL and 
h are the same as in (a)) with the effect of reionization taken into account 
(solid curve). The dashed curve represents the mass fraction of stars when 
the effect of reionization is not taken into account. 

(c) Evolution of the radiation energy emitted from stars (dashed curve) and 
from quasars (dashed-dotted curve). The solid curve represents the sum. 

(d) Comoving number density of quasars plotted against redshift for — 
lOVear and cqFq = 10"^ 

Fig.2 Critical temperatures for Tree = tc, ts — tc, Tree = H~^, tq — H~^, tc — 

plotted as a function of redshift. Region for Tree < "^c, or tq < tc, is 
needed for high energy photon emission from the collapsing objects. Another 
requirement that Tree < Tune, which is estimated from Fig. 6 below, is also 
shown. There are no regions that satisfy all the requirement. 

Fig. 3 (a) Evolution of the fraction of HI (sohd curve), HII (dashed curve) and Hc- 
parameter (dashed-dotted curve) shown for Q — 1, VLb = 0.05, h — 0.5 and 
/ = 0.02. 

(b) Evolution of metalicity Z as, a, function of redshift. 

(c) Evolution of the fraction of Hel (solid curve). Hell (dashed curve) and 
HeHI (dashed-dotted curve). 

(d) Evolution of electron temperature T'e(K), baryonic Jeans mass Mb,j, and 
the critical mass Mb^c above which the molecular cooling time is shorter than 
the dyanmical time. Mb,j and Mb,c are in units of Mq. 

Fig.4 Residual photon spectra at z = 5 (sohd curve), z = 3 (dashed curve) and 
^ = (dashed-dotted curve) for Q = 1, Qb = 0.05, h = 0.5 and / = 0.02. 



The dotted curve that closely follows the curve for z = 5 stands for the 
spectrum {z — 5) when the effect of absorption due to intergalactic H and He 
would be switched off. 

Fig. 5 Evolution of the ionizing flux from stars and quasars shown for fl — 1, fls — 
0.05, h — 0.5 and / = 0.02 (solid curves). Dashed curves represent the 
ionizing fluxes which are obtained when the effect of absorption due to H and 
He is switched off. 

Fig. 6 Condition for recombination cooling to dominate over line cooling for Q = 1, 
Qb — 0.05 and h — 0.5. The region below the sohd line shows the parameters 
for which recombination cooling dominates. The dashed curve shows the ratio 
n(HI)/n(HII) from ionization equilibrium. 
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